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ABSTRACT 

We use autocorrelation to investigate evolution in flow fields inferred by apply- 
ing Fourier Local Correlation Tracking (FLCT) to a sequence of high- resolution 
(0.3 "), high-cadence (~ 2 min) line-of-sight magnetograms of NOAA active re- 
gion (AR) 10930 recorded by the Narrowband Filter Imager (NFI) of the Solar 
Optical Telescope (SOT) aboard the Hinode satellite over 12-13 December 2006. 

To baseline the timescales of flow evolution, we also autocorrelated the magne- 
tograms, at several spatial binnings, to characterize the lifetimes of active region 
magnetic structures versus spatial scale. Autocorrelation of flow maps can be 
used to optimize tracking parameters, to understand tracking algorithms’ sus- 
ceptibility to noise, and to estimate flow lifetimes. Tracking parameters varied 
include: time interval At between magnetogram pairs tracked, spatial binning 
applied to the magnetograms, and windowing parameter a used in FLCT. Flow 
structures vary over a range of spatial and temporal scales (including unresolved 
scales), so tracked flows represent a local average of the flow over a particular 
range of space and time. We define flow lifetime to be the flow decorrelation time, 
r. For At > r, tracking results represent the average velocity over one or more 
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flow lifetimes. We analyze lifetimes of flow components, divergences, and curls 
as functions of magnetic field strength and spatial scale. We find a significant 
trend of increasing lifetimes of flow components, divergences, and curls with held 
strength, consistent with Lorentz forces partially governing hows in the active 
photosphere, as well as strong trends of increasing how lifetime and decreasing 
magnitudes with increases in both spatial scale and At. 


1. Introduction 


Estimates of photospheric velocities can be combined with the ideal Ohm’s l aw to dete r- 


mine the fluxes of magnetic energy and helicity across the photosphere (e.g., IChael 12001 


mine the nuxes ot magnetic energy and 
Demoulin and Bergeil2QQ.‘jl: Schuckl20odh 


These huxes probably play important roles in driv- 
ing the corona to hare (IWelsch et ah, 112009 1 and to launch coronal mass ejections (CMEs), as 
well as coronal heating fjTan et al.ll2007T) . Sever al techniques can be used to estimate p hoto- 


spheric h ows, including tra cking methods fe.g.. iFisher and Welsc. 


tracking fjPotts et all 120041) and feature tracking (IDeForest et al 


2008: ISchuckl 120061) . ball- 
2007 7 h . Techniques for es- 


timating photosph eric velocities, however, are imperfect (IRieutord et al.ll200ll; IWelsch et al. 
20 071; ISchu ck 2008]), so efforts to improve estimation techniques are ongoing. 


Autocorrelation of image sequences has been used to quantify evolutionary timescales of 
structures in the solar atmosphere. In this approach, a measure of correlation (e.g., linear 
or rank-order correlation) is computed between pixel values at some initial time T and the 
corresponding pixel values at a later time ti + AT, where we refer to AT as the lag time. 
As image features evolve with in creas ing; AT, measures of correlation computed this way 
typically drop. For instance, iTritschler et al.l (120071) investigated inter network and n e twork 


fine structure of the quiet Sun at disk center in the Ca II K line. Welsch et al. (2009) 
autocorrelated 96-minute-cadence, line-of-si ght (LOS) magne tograms of a few dozen active 
regions imaged by the MDI instrument (jScherrer et ah Il995h in full-disk mode (with ~ 2" 
pixels) over lag times ranging from the data’s nominal 96-minute cadence to several days. 
They found decorrelation times for the magnetic held, defined as a drop in the autocorrelation 
to 1/e ~ 0.37, on the order of 72 hours, with considerable variation between active regions. 


Welsch et al. (2009) also estimated horizontal hows in their active region sample by applying 
tracking methods to pairs of magnetograms separated in time by a tracking interval, At, 
equal to the data’s 96- minute cadence. The hows were independently est im ated by using 
two distinct methods, Fourier local correlation track ing (FLCT; iFisher and Welschl 120081) 
and the differential affine velocity estimator (DAVE; Schuck 20061 k Apodization windows 
of 8 and 9 pixels, respectively, were used, so motions on smaller scales were averaged away 
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over the apodization window. By autocorrelating the x- and ^-components of the estimated 
flows, they found flow lifetimes of two to three 96-minute steps on this ~ 12 Mm spatial 
scale. 


The autocorrelation approach employed by Welsch et al. (]2009i) was useful for estimating the 
lifetime of magnetic held structures and hows on the particular spatial scales they studied. 
But what are the lifetimes of magnetic and how structures on different spatial scales? How 
do how lifetimes vary as a function of length scale? And what are the lifetimes of how 
properties, such as vortical or converging motions? From the general properties of turbulent 
hows like those operating at the solar photosphere, we expect structure in how helds to be 
present over a range of spatial scales, with a range of lifetimes. 

Further, little effort has been undertaken to understand how variations in tracking cadence 


At, noise, and choice apodization window si ze affect how in 


these quantities influence estimated hows? Abramenko et al 


erence. 


How do variations in 


20111 used feature tracking to 


follow bright points over a range of At’s, and characterized their motion as super-diffusive^ 


The average speeds implied by their diffusive model decreased with increasing At. Chae et al. 


(2000) applied an LCT algorithm to chromospheric H -a images, and noted that At must be 
long enough for displacements to be large enough to be detected by LCT, but not so long 
that displacements are larger than the apodization window. They also found tha_t_ choosing: 


window sizes that were too small increased random velocities due to noise. Chae et a 


(1200 41) found slower average speeds for longer At’s. More recently, Verrna and Denkerl (20111 


conducted an investigation of the effects of parameter selection on hows inferred by applying 
local correlation tracking (LCT) to high-cadence, high-resolution (0.11 "pixel scale) G-band 
images of a held of view containing sunspot umbrae, penumb rae, an d quiet- Sun regions, 
observed by the broadband filter im ager (BFI) of the SOT (ITsuneta et al.ll2008l) aboard the 
Hinode satellite (IKosugi et al . 2007). The evolutionary time scales and spatial scales of the 
underlying image structures motivated their selection of tracking intervals between 15 and 
200 s, and apodization windows of 600 - 2400 km. They found decreases in how speeds both 
as the apodization window size increased, and as how maps were averaged over longer time 
periods. Because photospheric magnetic hux is long-lived, tracking methods can be applied 
to magnetograms separated in time by much larger tracking intervals — e.g., longer than an 
hour, as done by Welsch et al. 020091) . 


Here, we address these and other topics by using FLCT to track a nearly 13-hour sequence 
of high-resolution (0.16" pixels, with a ~ 0.32" diffraction limit), rapi d cadence (~ 2 min) 
LOS magnetograms of NOAA AR 10930, recorded by the SOT /NFI (ITsuneta et al. 20081: 
Suematsu et al.l 120081: iTchimoto et al] 120081 : IShimizu et al.l 120081) instrument aboard Hinode 
over 2006 December 12 - 13. This active region was the source of an X-class hare and 
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CME during our tracking interval: the flare be gan ar ound 13-Dec-2006 02:15UT and peaked 


around 02:40 UT in GOES 1-8 A X-rays, Tan et ah ( 2009 ) also investigated flows in this 


active region, focusing on the relationship of penumbra! flows and evolution to the flare, and 


Schriiver et al.l ( 20 080 conducted a comparative study of coronal magnetic field extrapola- 


tions derived from vector magnetograms of this active region. 


This paper’s primary goals are: first, to investigate the effects of choices of tracking parame- 
ters on inferred flows; and second, to characterize the lifetimes and other properties of flows 
as functions of spatial scale and magnetic field. We begin by describing the magnetogram 
dataset in more detail (§2]), then describe the methods we used to estimate flows (1J3]), present 
our analysis of the lifetimes of fields (1J4]) and flows ($5]), and conclude with a brief summary 
of our results and their implications (§HD- 


2. Hinode SOT/NFI Magnetograms 


We studied a sequence of Hinode SOT/NFI Fe I 6302 Ashuttered magnetograms of AR 
10930 with 0.16" pixels, created from the Stokes V/I ratio in Level 0 data, recorded between 
12-Dec-2006 at 14:00 and 13-Dec-2006 at 02:58. The USAF/NOAA Solar Region Summary 
issued at 24:00 UT on 12-Dec-2006 listed AR 10930 at S06W21, meaning it was near disk 
center during this interval. At this wavelength, the diffraction limit of the Solar Optical 
Telescope’s 50 cm aperture is 1.22A/a ~ 0.32". We show initial and final magnetograms 
from the sequence in Figure [IJ At that stage of the Hinode mission, a bubble present within 
the NFI instrument degraded image quality in part of the field of view. In these frames, 
the affected pixels are in an area between from about 300 - 900 in x and above about 900 
in y. To account for the bubble, we exclude pixels above y = 848 in all our analyses of 
flow properties and lifetimes (although this region is retained in some images shown in this 
paper). 


To convert the measured Stokes / and V signals in to pix el- averaged flux densities, we applied 
the approximate calibration used by Isobe et ah (120071 ) . 

C v 


Blos = 


x 10000 (Mxcm ) , 


( 1 ) 


0.798C7 

where £>los is the estimated LOS flux density, Cj and Cy are counts in the / and V images, 
respectively, and CV/0.798 gives the circular polarization. As noted by Isobe et al. (120071 ). 
this linear scaling is only approximate. In particular, in the dark cores of umbrae, the 
decrease in / introduces nonlinearity into the V/I ratio. This effect can produce spurious 
weak field regions in umbrae. A more complex, nonlinear calibration would be required to 
overcome this artifact. Since our primary goal is to investigate the time evolution of structure 
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in this magnetogram sequence, we have not pursued efforts to correct such artifacts in the 
V/I magnetograms. 

A detail ed spectroscopic ana lysis of fields in this acti ve region using Hinode SpectroPolarime- 
ter (SP nTsuneta et al.ll2008l) data was undertaken bv ISchriiver et al.l (120081) . who created one 
vector magnetogram with approximately 0.63 "pixels near the middle of our tracking interval, 
around 21:00 UT. Reprojection employed in their reduction procedure limits the ability to 
directly apply those results with our magnetogram data set, but, as described in Appendix 
we manually co-registered and resampled this SP vector magnetogram to map vertical 
magnetic field B z and field strength |B| from this magnetogram onto the (4 x 4)-binned NFI 
magnetogram closest to the midpoint of the 45-minute SP scan (at 21:52 UT). 

The AR was at about S06W21 toward the end of our observations, and NFI’s field of view 
(FOV) subtends about 10 heliocentric degrees at disk center. These facts imply cosines 
between the line of sight and the local vertical direction within the AR ranging from 0.88 to 
0.96, corresponding to variations in pixel scales from foreshortening of a few percent across 
the FOV. Over the earlier part of magnetogram sequence, upon which much of our analysis is 
focused, these variations are slightly smaller. Consequently, we take these modest distortions 
as acceptable, and opt not to reproject the magnetograms prior to tracking. Reprojection is, 
ho wever, certa inl y app ropriate for tracking studies over much larger holds of view, e.g., that 
of Welsch e t al. (120091) . who estimated velocities simultaneously across fields of view out to 
45° from disk center. 

Apart from three abnormally large time steps of approximately 600 s each, and two relatively 
small time steps of 26 s each, the 380 frames we analyze have a cadence of 121.4 ± 1.2 s. 

To study frame-to-frame changes in the magnetic held, we must hrst remove artiheial evolu- 
tion due to whole-image shifts (rigid motion of each image) between frames, due primarily 
to variations in spacecraft pointing. To accomplish this, we hrst determined the frame-to- 
frame shift between the central regions of each pair of images (cropping the outer 10% of 
pixels along each edge), using FFTs to compute the frame-to-frame cross correlation func- 
tion and a second-order Taylor series approximation to es timate the displacement vector of 
the function’s peak, to sub-pixel accuracy (iFisher and Welschl 120081) . We then used Fourier 
interpolation to estimate pixel values implied by a shift of the negative of the cumulative 
displacement from the initial frame in the sequence. Any wrapped pixels arising from this 
shift, due to the assumed periodicity, were then zeroed (as may be seen in a thin strip along 
the bottom of the right panel of Figured]). In Figured] we show the cumulative shifts in 
x and y over the duration of the magnetogram sequence. Large shifts in image alignment 
occur with a n ap proximate periodicity similar to that of Hinode’s 98-minute orbital period. 
(iKatsukawa et alJ[2010 ). 
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Fig. 1. — Left: Initial Hinode/NFI LOS magnetogram (white is positive flux, black is negative) from the 
13-hour data sequence we investigated, at full resolution (0.16 "pixels) with the saturation level set at ±500 
Mx cm~ 2 . Right: Final magnetogram in the sequence. 


In Figure [21 we also show the total unsigned flux and negative of the net flux (there is more 
negative LOS flux than positive; but the right plot axis is positive, so we changed the sign 
of the net flux), which both exhibit periodicities on a timescale similar to that of the whole- 
frame shifts. Analysis of Fourier spectra of the time series in pixels (not shown) shows a 
significant excess power near the orbital frequency, but little evidence of helioseismic P-mode 
leakage into the magnetogram signal. Given the large Doppler shift from orbital motion, 
observed periodicities could arise from leakage of Doppler shifts induced by spacecraft orbital 
motion into the Stokes V measurement. Thermal variations in the instrument might also be 
responsible. 


As discussed in more detail below, local fluctuations in held strength from noise can mimic 
changes from true magnetic evolution, implying estimated hows will be inaccurate in regions 
where the measured helds are n ear t he magne tograp h noise level. We estimated this noise 
level using the procedure of IHagenaar et ah (119991 ) . First, we spatially binned the data 
2 x 2 to 0.32"-wide pixels, to approximate the SOT diffraction limit of 0.3". Next, we 
created histograms of pixel values in each frame in the sequence, in bins of 3 Mx cm -2 . 
Then, assuming the core of the distribution in signals arises purely from noise, we ht a 
Gaussian, exp[— (B — B 0 ) 2 /(2a 2 ), to the central region of the histogram, ±10 Mx cm -2 . 
Over the course of the tracking interval, htted values of the width parameter a for single 
magnetograms ranged from a high near 17 Mx cm~ 2 at the start of the run to lows near 
7 Mx cnr 2 toward the end of the run. For simplicity, we therefore take 15 Mx cm -2 as a 
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Fig. 2. — Left axis: Cumulative frame- to- frame shifts (in pixels) in the x (red) and y (blue) directions be- 
tween pairs of magnetograms in the NFI sequence that we analyzed, with respect to the initial magnetogram. 
Right axis: total unsigned flux (dashed) and the negative of net flux (solid). Variations in flux probably 
arise from crosstalk due to Doppler shifts induced by spacecraft orbital motion into the Stokes V signal. 

uniform estimate of the noise level. 


3. Estimating Flows 

Assuming magnetic diffusi vity is negligibly small over the time interval At between a pair of 
successive magnetograms, Demou lin a n d Ber ger (2003) argued that the apparent motion of 
magnetic flux from one magnetogram to the next over At could be described by a “footpoint 
velocity” u, that is related to the plasma velocity v by 


u B. n = vB, - v„B h 


( 2 ) 


where n denotes the unit vector normal to the surface, the n subscript denotes the normal 
component of a vector, and the h subscript denotes the horizontal components of a vector, 
which are perpendicular to n. This assumption can be used to express the normal component 


of the magnetic induction equation in terms of a continuity- like, finite-difference equation, 


A B r 
At 


+ V/i • (u B n ) — 0 , 


( 3 ) 


where A B n = B n (tf ) — B n (ti ) is the change in the normal magnetic held between the initial 
time t t and the final time tf, and V is approximated with finite differences in the spatial 
coordinates orthogonal to the normal direction. The analogy with the continuity equation 
is not precise, because unlike number or mass densities, which are always positive, magnetic 
flux can be positive or negative. So, unlike total particle number or total mass in a volume, 
neither total signed nor total unsigned flux is conserved in a magnetogram: existing opposite 
polarity fluxes can cancel, or new bipolar flux can emerge. 

Dcmoulin and Berger ( 20031) argued t hat tracking metho ds ap plied to magnetograms esti- 
mate u, instead of w h (cf., Chad 1200111 . but S chuck (120051 1200811 argued that some tracking 
methods return a biased estimate, u, with the contribution from weighted more heavily 
than in the true u. The term “footpoint velocity” implicitly refers to extensions of magnetic 
flux above the magnetogram surface. Welschl (2006) used the alternative “flux transport 
velocity” to refer to u without reference to held structure above or below the magnetogram. 
For purposes of tracking, we assume the LOS held, Be, derived from the NFI observations, 
approximates the normal held B n . 


We used the FLCT code (jWelsch et al.ll2004J: iFisher and Welschll2008l) to estimate hux trans- 
port velocities for this dataset. The FLCT code uses Fourier cross-correlation of B n (ti ) with 
B n (tf ) near a given pixel, say one centered at ( Xj,fjk ), to determine the two-component dis- 
placement AR of magnetic held structure between the magnetograms near {xj, y fc ); division 
by At then results in a velocity estimate at ( Xj,yk )• The correlation is localized by hrst 
weighting each magnetogram by a Gaussian windowing function, exp(— r“j k /a 2 ), where is 

distance from ( Xj,yk ), and a is a free parameter. Optionally, low-pass filtering, with a Gaus- 
sian roll-off, can be applied to the Fourier spectra prio r to performing; the cross-correlation. 
In tests applying shifts to quiet-Sun magnetograms, Fisher and Welsch (120081) found roll-off 
values of 0.2fc max - 0.5/c max improved recovery of the applied shifts, where fc max is the highest 
wavenumber each dimension. In this work, we used a roll-off wavenumber of 0.5 fc max , which 
corresponds to relatively weak Fourier filtering. 

Like FLCT, other tracking methods typically apply an apodization or windowing function, 
the spatial extent of which is a free parameter, to the initial and final magne tograms, and 
find a velocity that extremizes some functional over this window ( iSchuckl 2006). Given data 
with sufficient cadence, the interval At between pairs of magnetograms to be tracked can also 
be varied. Tracking methods therefore estimate velocities from the effects of actual velocities 
averaged both spatially over the windowed region and temporally over the time interval At 
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between magnetograms. This averaging implies that flows on spatial scales smaller than the 
windowing function, and with lifetimes shorter than At, cannot be recovered in the flow 
estimation process. 

Beyond these algorithmic aspects of the tracking process, it can be also seen that noise 
and other artifacts in magnetogram measurements will contribute to the difference A B n in 
equation (J3]), implying such effects can influence estimated velocities. 

These considerations motivate efforts to understand the interplay between windowing and 
tracking cadence and flow spatial scale and lifetime, as well as the effects of magnetogram 
noise. To investigate these relationships, we estimated flux transport velocities from pairs of 
magnetograms that were separated in time by eight varying intervals, with At G {2, 4, 8, 16, 
32, 64, 128, 256} minutes, that were both unbinned and binned to six varying macropixel 
sizes Ax G {2, 4, 8, 16, 32, 64 } in terms of the data’s original 0.16" pixels, with four 
choices of apodization parameter a G {2, 4, 8, 16 }. Pixel values were simply averaged 
when rebinned, which does not exactly mimic how magnetographs with either worse spatial 
resolution or poor weak-field sensitivity would detect magnetic flux. 

Since we are interested in characterizing the effects of noise on estimated flows, we tracked 
all pixels above a threshold of 15 Mx cm -2 , which is both near our estimate of the noise 
level, and includes the spurious weak-field regions in the sunspot cores. 


4. Lifetimes of Magnetic Structures 

To provide context for the time evolution of flows on varying spatial scales, we autocorrclated 
the magnetograms themselves, at several spatial resolutions: with 0.16" pixels, and then 
with data rebinned to macropixels Ax G (2, 4, 8, 16, 32, 64} pixels on a side. While each 
magnetogram is a 2D array, we computed the correlations as if both were ID arrays. All 
pixels outside of the bubble region were used, not just those above the noise threshold used 
for tracking. The colored curves in Figure [3] show rank-order correlation coefficients for 
the autocorrelations at each spatial resolution for increasing lag time. In addition, the black 
dashed and solid curves show frame-to- frame linear and rank-order correlations, respectively, 
for the 2x2 binned data. Note the limited range of the y- axis. 

The high frame-to-frame correlations imply that the magnetic field does not evolve much 
from one magnetogram to the next. It can also be seen that frame-to-frame linear correla- 
tions hover around unity, but frame-to-frame rank-order correlations are significantly lower, 
implying rank-order correlation is a more sensitive test of differences between magnetograms. 
This plot also shows that the lifetime of this active region’s magnetic structure, over the en- 
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B LO s Autocorrelation vs. Bin Size 



Fig. 3. — Colored curves show rank-order autocorrelation coefficients with increasing lag times, for magne- 
tograms at several spatial binnings: red is for 0.16 "-pixel data; orange, yellow, green, aqua, blue, and purple 
are for data rebinned 2 x 2, 4 x 4, 8 x 8, 16 x 16, 32 x 32, and 64 x 64, respectively. The black dashed and 
solid lines show frame-to- frame linear and rank-order correlation coefficients, respectively, for data binned 
2x2. Note the limited range of the vertical axis. Downward spikes before 500 minutes arise from frames 
with many cosmic ray hits. The effects of the X-class flare can be seen as small blips near 743 minutes. 

tire magnetogram field of view, is much longer than the duration of the dataset that we 
analyze. 

Several artifacts are also present in these time series. Downward spikes in correlations 
correspond to frames in which isolated pixels contained large fluctuations, probably due to 
hits by cosmic rays. Despiking the magnetograms can remove these spikes, but here we are 
interested in the effects of such artifacts on autocorrelations. Larger binnings average out 
these effects. In terms of lag times plotted here, the start time of X-class flare that began 
around 02:28UT on 2006-Dec-13 corresponds to 743 minutes, and its effects can be seen as 
small wiggles in some autocorrelation curves near the right end of the plot. 

Distinct types of magnetic regions — notably, sunspots, plage, and quiet Sun — all influence 
the autocorrelations in Figure |3] By investigating the lifetimes of magnetic structures on 
smaller spatial scales, we can better understand the lifetimes of magnetic structures in each 
type of region. Accordingly, in Figure HI we show autocorrelations (red curves) for -Blos 
in (32 x 32)-binned subregions of (2 x 2)-binned data (so each subregion is (64 x 64) full- 
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resolution pixels). In the same figure, we also show one-parameter fits (blue curves) to the 
decorrelation in each subregion, assuming the autocorrelation function decays exponentially 
with lag time AT, exp(— A T/r), with the decay constant r as the only free parameter. In 
each subregion, we take the fit decay constant, r, to be the lifetime of the structure in that 
subregion. We define the occupancy of a subregion to be the percentage of pixels above our 
15 G noise threshold. Only pixels with a median occupancy, over the dataset’s duration, of 
at least 20% were fit. Background gray contours in the figure show 50 G and 200 G levels 
of | B Tj os | from the initial magnetogram (in full- resolution pixels), so variations in lifetimes 
can be understood in the context of active region structure. 


We then investigated the dependence of magnetic lifetime on field strength. To start, we 
computed a measure of subregion field strength by boxcar averaging flux densities from 
|i?Los| from the five initial magnetograms in our sequence at (2 x 2) resolution, then taking 
the absolute value, then computing the spatial average of this unsigned, boxcar-averaged 
flux density over each subregion. We chose to use this measure of initial field strength to 
address the question: “For a given current field strength, how long will a structure persist?” 
The relatively short fitted lifetimes of some magnetic regions suggests that averaging over 
longer intervals — e.g., the length of the entire run — would imply averaging over several 
magnetic turnover times. It can be seen in Figure [4] that some of the fits represent the 
evolution poorly, but many fits match the observations well, given the single free parameter. 
Of the 195 subregions with at least 20% above-threshold occupancy, only four have fitted 
lifetimes that are less than the two- minute nominal cadence (see, e.g., the two subregions 
to the left of the top-right subregion). We ignore these too-short fits as pathological. In 
Figure [5], we plot the fitted lifetimes as a function of these spatially and temporally averaged 
subregion fields (blue +’s). A linear fit between the logarithms of each returns a power- 
law exponent of 0.66. The solid blue line shows the fit, while the dashed blue lines reflect 
uncertainty in the fitted slope, based upon y 2 computed from the fit. This trend is consistent 
with the idea that convection disr upts magnetic field structures, but stronger magnetic field s 
more strongly inhibit convection fjTitle et al.l Il992l : iBerger et al.l Il998l: IWelsch et al.l 120091 1. 
so structures persist longer in stronger-field regions. Despite the trend of increasing lifetime 
with field strength, we note some high-field-strength regions have short lifetimes. Relatively 
rapid evolution in these high-field strength regions might play some role in flare activity, but 
pursuing this question is beyond the scope of this paper. 


To investigate the dependence of this field-strength-versus-lifetime relationship on spatial 
scale, we repeated this procedure with the subregion size reduced to 16 2 (2 x 2)-binned 
pixels (autocorrelations and fits at this resolution are not shown). Lifetimes as a function of 
subregion- averaged field strength are plotted in red in Figured excluding 28 subregions with 
fitted lifetimes < 2 min (of 733 subregions with at least 20% above-threshold occupancy). 
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A linear fit between the logarithms of held strengths and lifetimes returns a power-law 
exponent of 0.69 in this case. The solid red line shows the fit, while the dashed lines reflect 
uncertainty in the fitted slope. The trend of lifetime versus held strength did not change 
much for subregions smaller by a factor of two. 


Recalling that our estimates of |Rlos| suffer from non-monotonicity in the Stokes V/I ratio, 
we also characterized the lifetimes of NFI structures as functions of \B Z \ and |B| from the SP 
data, in subregions averaged to the same spatial resolution (16 2 regions of 2 x 2 pixels). Note, 
however, that the SP data were recorded near the middle of our observing run (between 20:30 
and 21:14 UT), while the magnetic lifetimes were computed from the beginning of our run. 
Also, our treatment takes no account of spatial filling factors of magnetic structures within 
resolution elements. Despite these uncertainties in the SP held strengths, hts of lifetimes as 
a function of \B Z \ (not shown) were similar to those for |Slos|- The solid green line in Figure 
[5] shows the htted slope (dashed lines rehect uncertainties in the htted slope) of lifetimes 
as a function of |B|. The slightly weaker dependence of lifetime on |B| compared to |I?los| 
might arise from either: vertical magnetic helds inhibiting convection more efficiently than 
horizontal magnetic helds; or rapid magnetic evolution near polarity inversion lines (PILs) 
of the normal held, regions where the magnetic held is predominantly horizontal and where 
hux emerges or submerges; or both. We note that the range of | B | is greater than for either 
of the other two qu ant iti es, and t his longer run will naturally decrease the htted slope for 


the same rise. IDemoulin and Bergen ( 20031 ) suggested that in regions with predominantly 


horizontal magnetic helds, apparent motions of vertical magnetic hux could be rapid. 


We remark that autocorrelation analyses from observed magnetograms like those presented 
here could be compared to autocorrelati on s fro m synt hetic magnetograms extracted from 


simulations of magnetoconvection (e.g 
between the observations and simulations 


Rempel 2011), as a test of statistical consistency 


5. Flow Lifetimes 

In this section, we hrst consider selection of tracking parameters, and then investigate life- 
times of estimated hows versus spatial scale and magnetic held strength. 


5.1. Selection of Tracking Parameters 


Welsch et al. 02 004) n oted that hows derived by tracking alone only approximately solve 


equation (J3]). Schuck 020061 ) suggested that, in the presence of magnetogram noise, exactly 
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solving the induction equation is probably not optimal. W elsc h et al. (1200711 suggested that 
for fixed At (i.e., imposed by available data), cr could be chosen to optimize agreement 
between the observed AB n /At and the estimated • (uf? n ). When observations are fre- 
quent enough that At can also be varied, considerations beyond optimizing consistency with 
equation ([3]) are useful. 


In addition, Welsch et al. (2 011( 1 demonstrated that enforcing strict consistency with equa- 
tion ((2]) in deriving flows is probably unwise. Using 96-minute cadence, full-disk magne- 
tograrn sequences observed by MDI dScherrer et ah 1995 1. they characterized magnetic evo- 
lution in three active regions near disk center. In addition to using FLCT to estimate flows, 
they also numerically derived ideal electric fields to exactly solve Faraday’s law, from a 
Poisson equation, with a source given by their estimate of the radial magnetic held evo- 
lution \ABr/ A t], for the inductive part of the electric held (IWelsch et al.l 120071 ) . Where 
Welsch et ah (2009) found how helds derived by either FLCT or DAVE from 96-minute ca- 
dence magnetogr ams to be c orrelated from one 96-minute interval to the next, the “Poisson 
how holds” that Welsch et al. (120111) derived to be exactly consistent with ABr/ At decorre- 
lated from one frame to the next. Assuming random huctuations (due to, e.g., noise present 
in the estimated hows), it can be shown that, if evolution in the how is negligible (so the 
only difference between the how maps is due to random errors in the estimated hows), then 
even a relatively high correlation coefficient (e.g., 0.9) implies significant huctuations present 
(respectively, 33% of the variance in the how component). It can also be shown that even 
allowing for evolution in the how only moderately reduces the fraction of decorrelation that 
must be ascribed to random huctuations in the how estimates. The fact that FLCT and 
DAVE hows did not decorrelate between frames, but the Poisson hows did, implies the latter 
contained significant random huctuations, and were therefore probably strongly influenced 
by magnetogram noise. This suggests that using equation (J3]) as the sole guide to determining 
tracking parameters has limitations. 


We therefore consider the effect of the time interval At between each pair of tracked mag- 
netograms on the relationship between hux transport velocities u from tracking and plasma 
velocities v in equation (J2]). Essentially, if At is either too short or too long, then equation 
((2]) is invalid. To explain why, we consider each case separately. 

When the time difference At between magnetograms is too short, very little magnetic evolu- 
tion has occurred (as in the black curves in Fig. [3]), so most of the difference AB n between 
initial and final magnetograms that appears in equation ([3]) is due only to noise. Conse- 
quently, hux transport velocities inferred from AB n / At will be spurious. For instance, a 
transient spike in the held strength of a single pixel in a weak unipolar region would lead 
to inference of a high-speed converging how surrounding that pixel, consistent with the sud- 
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den field increase, followed by a diverging flow, consistent with the field decrease after the 
spike disappears. In general, these noise-induced flows will be unrelated to actual plasma 
velocities. Consequently, if At is too short, then noise is a large component of magnetogram 
difference A B n , and we refer to flows estimated in this regime as noise-dominated. 

When the time difference At between magnetograms is too long, displacements of magnetic 
flux might exceed a pixel length, meaning the finite difference approximation used in equation 
© could be inaccurate. In this regime, which we refer to as displacement-dominated, flux 
transport velocities estimated by tracking are more properly described as time averages of 
true flux transport velocities, which acted and evolved over timescales shorter than At. It 
should be noted that, although estimated flows in the displacement-dominated regime do 
not correspond directly to plasma velocities, they still have a clear physical meaning: they 
represent the time-averaged flow over timescales of a flow lifetime or longer. We expect that 
motions derived for At longer than a flow lifetime will underestimate true flow speeds: the 
true path of a feature executing zig-zag motion can be resolved with short enough At, but 
if At is too long to resolve each step on the path, only the (shorter) net displacement is 
recovered (cf., the full path length traversed over that time interval). Hence, a given At can 
average over the effects of shorter-lived flows. 

It would be convenient if an ideal choice of At could be made a priori. As previously noted, 
however, flows with a range of speeds, and time and length scales are present at the solar 
photosphere (e.g., granular flows with speeds ~ 1 km sec -1 with lifetimes of a few minutes 
and ~ 1 Mm length scales, and supergranular flows with speeds a few times 10 2 m sec -1 with 
lifetimes of about a day, and length scales of ~ 10 Mm). Poor choice of At can, however, 
be determined by post facto analysis of tracking results, meaning trial runs can be used to 
identify suitable timescales. 

We first consider the noise-dominated regime, in which correlations between each velocity 
component from successive flow estimates will be small. Frame-to-frame correlations can 
be improved by increasing At — t f — U, and by averaging successive magnetograms around 
ti and tf to synthesize initial and final images. Figure [6] shows frame-to-frame, rank-order 
correlations for the x- and y-components of estimated velocities, with a fixed at 4 pixels but 
for several Af’s, using either unaveraged (dashed) and averaged (solid) magnetograms for 
B(ti) and B(tf). Low correlations for short values of At imply that noise has degraded the 
velocity estimates. For At = 2 min, the flow fields are even weakly anticorrelated. Increas- 
ing At clearly increases frame-to-frame correlations, because real magnetic field evolution 
becomes more significant than spurious evolution from noise. Saturation occurs, however: 
doubling At from 8 to 16 minutes did not greatly improve frame-to-frame correlations. As At 
increases further, into the displacement-dominated regime, frame-to-frame correlations can 
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approach unity, since displacements averaged over very long time intervals do not change 
much if the intervals’ endpoints are shifted slightly. To average the magnetograms, a 3- 
magnetogram boxcar was used for At = 2 and 4 min, while a 5-magnetogram boxcar was 
used for At = 8 and 16 min. For cadences of 2, 4, and 8 minutes, there is overlap in data 
used for B(U ) and B(tf). This probably leads to blurring of features from their motion y>ver 
the s umming interval, and could also lead to correlated noise in B(ti) and B(tf ) (Schuck 
20061 ). The improvement in frame-to- frame correlations from averaging magnetograms is 


substantial. Linear correlation coefficients (not shown) were both lower and less consistent 
across frames. We therefore employ rank-order correlations to compute all subsequent auto- 
correlations. We provide concrete examples of these effects in Figures [7] and |8j In the former, 
successive flow fields computed with At = 2 min show very little consistency between flow 
fields. In the latter, successive flows computed with At = 8 min are much more consistent. 
Generally, frame-to-frame correlations are higher with larger values of a (not shown). This is 
consistent with FLCT detecting flows on larger spatial scales with larger cr, and larger-scale 
flows persisting longer in time. We discuss this issue in greater detail below. 


In the displacement-dominated regime, correlations between each velocity component from 
flow fields separated by lag times AT much shorter than At can be high, but such correlations 
will be very weak when AT is on the order of At or longer. This implies substantial flow 
evolution has occurred over the At between magnetograms. Phrased another way, At is 
longer than the estimated flows’ lifetimes. For estimated flows to reflect actual plasma 
velocities, however, At cannot exceed the flow lifetime. Displacement-dominated flows are 
still physically meaningful, however: they quantify the average effects of flows over the course 
of more than one flow lifetime. 


Having demonstrated the advantages of averaging, we use flows derived from magnetograms 
averaged with a five-step boxcar in our analyses of flow lifetimes below. For convenience, 
we will sometimes refer to flows inferred by FLCT as velocities, but the reader should keep 
in mind that FLCT’s flows should not necessarily be interpreted as estimates of plasma 
velocities at the center of the time interval At. 


5.2. Estimating Flow Lifetimes 

We define the lifetime r of a flow as the lag time AT at which a one-parameter fit to the 
autocorrelation drops to 1/e ~ 0.37. In Figure [9] we plot autocorrelations for flow estimates 
made with two choices of cr and several choices of At. For At = 64 minutes, the lifetimes of 
the x- and y- flow components are around 35 minutes for cr = 8, implying that flow estimates 
made with At = 64 minutes are displacement-dominated. In contrast, the lifetimes of cr = 8 
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flows estimated with At = 4, 8, and 16 minutes are all at least as long as At. Lifetimes of 
flows estimated with cr = 4 (dashed lines) are about half as long as a = 8 estimates for At = 
8, 16, and 64 minutes. 

The autocorrelations in Figure [9] also suggest that flows estimated on larger spatial scales 
- using a larger value of cr — have longer lifetimes. In fact, this is typical. In addition, 
longer flow lifetimes are also found when data are rebinned into larger macropixels, though 
there is a point at which rebinning does not increase decorrelation times. The top panel of 
Figure [TU] demonstrates that choosing larger values of cr or rebinning into larger macropixels 
are practically equivalent. The bottom panel of Figure [TO] implies that, when flows are 
displacement-dominated (i.e., they decorrelate at lag times AT < At), longer-lived flows 
can be identified in rebinned data. We note that increasing the windowing parameter cr can 
increase compute time for FLCT (and probably other tracking methods), while rebinning 
(holding a fixed) can decrease compute time, since fewer velocity estimates must be made. 

In Figure [IT] we show autocorrelations of u x and u y in (32 x 32)-pixel subregions of mag- 
netograms with 0.32 "pixels (so binned 2x2). These flows were derived with At = 64 
minutes and a = 4 pixels. Background gray contours in the figure show 50 G and 200 G 
levels of |-Blos| from the initial magnetogram (in full-resolution pixels), so variations in flow 
autocorrelations can be understood in the context of active region structure. 

In weak-held regions, few pixels are above the threshold for tracking, so very few velocities 
are estimated. Consequently, the autocorrelations are essentially random. We therefore only 
fitted the autocorrelations in high-occupancy subregions, which we define to be those with a 
median (over the dataset’s duration) occupancy of pixels above the tracking threshold of at 
least 20%, which gives large enough samples that the autocorrelations are stable. The fits to 
the magnetic held autocorrelation curves in Figure 0] convey how well the single-parameter 
exponential fits match the autocorrelations for a range of lifetimes and autocorrelation pro- 
files. 

In Figure fl2l we show fitted how lifetimes ( u x = red +’s and u y = blue x’s) for high- 
occupancy subregions of Figure dH as a function of the spatial average of the unsigned, 
boxcar-averaged LOS magnetic hux density in each subregion. Red A’s and green V’s are 
lifetimes of u x components of the how as functions of spatial averages of \B Z \ and |B|, respec- 
tively, and red dashed and green solid lines are hts to each. Most how lifetimes are between 
about 7 and 30 minutes, conhrming the expectation based upon the red, dashed curves in 
Figure El that these hows are displacement-dominated. There are statistically significant cor- 
relations both between how lifetime and held strength, and between how lifetime and median 
occupancy of above-threshold pixels. The rank-order correlation of held strength |B| with 
u x (and u y ) how lifetime, however, is a statistically significant 0.52 (and 0.51) for the 35 of 
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137 non- bubble subregions of SP data with occupancies of 1000 or more (of a possible 1024), 
implying the lifetime versus field-strength trend is independent of any occupancy-lifetime 
correlation. Linear fits to logarithms of lifetime as a function of held strengths return weak 
power-law exponents, below ~ 0.2 for each how component. For hows with At = 8 min- 
utes and a =4 pixels, correlations of lifetime with held strength are quite weak. For hows 
estimated with At = 256 minutes and a =16 pixels, htted power-law exponents between life- 
time and held strengths {I-Blos^ \B z \, | B | } range from 0.2 0.3, but held strength-lifetime 

correlations are weaker for u x than for u y . The east-west orientation of the PIL, combined 


with apparent rapid motion of vertical magnetic hux along the PIL flDemoulin and Berger 
20031 1. could produce more rapid evolution of east- west how components than north-south 
how components. 

Correlations between how lifet im e in a subregion and magnetic held strength in that subre- 
gion imply that Lorentz forces drive a component of photospheric evolution. To the extent 
that convection is the primary driver of photospheric evolution, and magnetic helds inhibit 
convection, there is no obvious reason why how lifetimes should be longer in stronger-held 
regions. But because magnetic structures persist for longer times cales than hows, it is likely 
that forc es acting on those he l ds, su ch as magnetic buoyancy (IParkerl 119551) or magnetic 
torques flLongcope and Welsch 20Q0) are also longer-lived. Consequently, it is physically 
reasonable that in stronger-held regions, magnetic evolution is partially driven by longer- 
acting forces than the pressure gradients arising from convection. 

How does how lifetime r depend upon how spatial scale? A crude estimate of the turnover 
time of a convective eddy — expected to be related to how life time — is the eddy length scale 


L, divided by a characteristic speed uq. As noted bv lBerger et al.l f )l998i ). however, “the term 
‘characteristic speed’ must be used only in relation to a given spatiotemporal resolution.” So 
as a preliminary step in characterizing how lifetime versus spatial scale, we hrst investigate 
how speed versus spatial scale. What do our results reveal about typical how speeds over 
varying spatial scales? We dehne a how’s spatial scale as the product of macropixel bin 
size and the windowing parameter a. As shown in Figure [TO], changes in both parameters 
seem to have similar effects upon how lifetimes. Because we estimated hows for several 
macropixel sizes and choices of cr, we therefore have multiple sequences of how estimates 
at a few spatial resolutions for each cadence At. We computed the average how speed at 
each scale over all pixels and time steps for which velocities were estimated. In Figure [GU 
we show average velocity as a function of spatial scale; the plot symbols are log 2 (bin size), 
in units of the data’s native 0.16 "pixels, to better illustrate changes in average velocity as 
a function of bin size at a given spatial scale. (Values of curls and divergences of the how 
pattern are also shown in this hgure, and are discussed below.) One basic pro perty is evident: 
average speeds are slower for longer tracking intervals At, as was noted by Chae et _ah 020041 ) 
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and V erma and Denke rl (2 0 111 ) for both relatively long and short At, respectively. (This is 
discussed further below.) For At shorter than 16 min, a verage speeds peak at the small end of 
the range of spatial scale s, as found in previous work (IBerger et al.l Il998c IChae et adj [2004J; 
Verma and Denkerl 2011). For longer At, however, average speeds peak at higher spatial 
scales: speeds at small-scales for large At are low. The large At implies that small-scale 
structures at U and tf have evolved substantially. While it is plausible that FLCT would 
infer random displacements in this case, and that averaging over these random displacements 
(all positive) would produce a large mean speed, we do not observe this. Apparently, for a 
given At, the FLCT code is primarily sensitive to flows on a particular spatial scale. Also, for 
shorter At, average speed at a given spatial resolution depends strongly upon spatial binning, 
while for longer At, average speeds seem insensitive to spatial binning. This probably reflects 
different effects from the combination of actual evolution and noise on the apparent magnetic 
evolution convolved with spatial binning and apodization by a. For instance, large values 
of a presumably average over the effects of actual evolution on sub-cx scales. If so, speeds 
should be higher with smaller a at a given spatial resolution, as we find. We remark that the 
combined effects of cadence with binning and apodization preclude straightf orward t esting 
our FLCT velocities for consistency with the Kolmogorov scaling (e.g., Thompson 2006 , 

p.68). 

We note that tests of FLCT with synthetic magnetograms from MHD simulations, in which 
the true vel ocities were known, showed that FLCT systematically underestimates speeds, 
by 25% flWelsch et al.l 120071 ) . so the absolute values of the average speeds we find here 
are probably systematically lower than estimates that would be found using other tracking 
codes. Further, because our velocities are derived from magnetic evolution, and magnetic 
flux both tends to inhibit convection and to concentrate in stable downflow regions between 
convective cells, it is reasonable to expect our average speeds to be systematically lower than 
tracking results derived from intensity images. We note that for At of 2 or 4 min, our average 
speeds at the smallest spatial scale are close to those reported by IVerma and Denkerl (2011) 
their Table 3). The speeds of FLCT flows derived by Welsch et al. ( 20091 ). with At of 


see 


96 min and a spatial scale of ~12 Mm were predominantly in the 50 - 100 m sec 1 range 
(see their Fig. 6), comparable to values we find for similar At and spatial scales here. 


In Figure [TT1 we show the average speed versus tracking cadence At, averaged over both 
space and time in all flow maps at a given spatial resolution and all spatial resolutions 
at each cadence. Effectively, this amounts to computing the average of speeds at each 
At in Figure [H2 — fitting the speeds in each panel to a horizontal line. Error bars in 
Figure [TTl reflect the standard deviation in speeds in each panel of Figure [TSJ A linear fit 
to the lo garit hms of cade nces and average speeds returns an power-law exponent of -0.34. 
Verma and Denker (20111) computed average speeds for a range of much shorter Af’s, and 
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over the eight-fold increase in At from 60 to 480 sec, their average speed fell by about a 
factor of two. This implies a scaling near -1/3, similar to ours. (But they found an increase 
in average speed between Af’s of 15 and 60 sec, and similar average speeds at 60 and 90 
sec.) They suggest decreasing speeds with increasing At might arise from the short lifetimes 
of features in the images they track. However, we find the same effect when tracking much 
longer-lived magnetogram structures, implying the limited lifetimes of features is probably 
n ot so lel y respon sible for decreasin g speeds with increasing At. We also note that both 
Hagenaar and Shine (2005) and Verma et al. ( 2011 ) found that individual, successive moving 
magnetic features (MMFs) followed preferential paths as they moved outward from sunspots, 
suggesting short-lived features can reveal longer- lived flow patterns. As noted above, a longer 
At can average over the effects of short-lived flows: displacements inferred with a longer At 
might not recover zig-zag motions that can be resolved with shorter At, with the shorter 
apparent path length at longer At implying a lower estimate of speed. 


Verma and Denker (20111) also averaged higher-cadence flow maps over increasing time pe- 
riods — 1, 2, 4, 8, and 16 hr — and found statistics of averaged-flow speeds (mean, median, 
10th percentile, maximum) all systematically decreased. Such averaging is not the same as 
tracking over longer At, as we have done, but the trends are qualitatively similar: a log-log 
plot of the averaging peri ods and mean speed from tracking granulation given in Table 1 of 
Verma and Denk er 120111) also shows a power-law decrease of mean speed with time between 
1-8 hr, though a linear fit to the logarithms returns an exponent near -0.18, about half 
our result. It is plausible that short-lived flows, resolved with shorter At, would primar- 
ily reflect the essentially random effects of granulation, and therefore average incoherently 
over long times. H owever , comparison of our large- At flow speeds with speeds inferred by 
Verma and Denker 1120111) with a shorter tracking interval At' and then averaged over longer 
time intervals similar to ours implies that shorter-timescale flows are averaging coherently. 
It is possible that the motion tracers — either the intensity pattern, or magnetic flux — used 
in one approach or the other (or eve n both!) do not faithfully represent the true plasma 
velocity. In this vein, Rie uto rd et al. 1120011 1 has argued that the motion of granules does 
not accurately reflect underlying velocities on spatial and temporal scales that are too short. 
This suggests flows inferred from magnetic tracking by flow maps derived with short tracking 
intervals but averaged over long At and flow maps derived with a tracking interval of At 
should be compared in a future study. 


We now return to the question of how flow lifetime r varies with spatial scale. As above, we 
define the spatial scale of an estimated flow as the product of macropixel size and windowing 
parameter a, and therefore have multiple sequences of flow estimates at a few spatial reso- 
lutions for each cadence At. For each sequence of flow estimates, we can estimate the flow 
lifetime r over the whole active region from a one-parameter fit to the autocorrelation as 


20 


function of lag time AT, assuming exponential decay, exp (— AT/r). In Figure |T5l we plot 
flow lifetimes estimated this way, as a function of spatial scale, for several tracking cadences 
At. In the shorter-scale region of each log-log plot, flow lifetime seems to scale, roughly, in a 
straight line with spatial scale. To characterize this dependence, we fit power-law exponents 
for each cadence over the range of seemingly straight-line log-log scaling (shown with a solid 
line on each plot). For longer cadences, flow lifetime appears to scale nearly linearly with 
scale length, while the dependence of flow lifetime on length scale seems weaker for shorter 
cadences. At larger length scales in each plot, the relationship between lifetim e and length 
scale becomes less coherent as binning and apodization are varied. Welsch et ah (120091 ) 
found lifetimes at ~ 12 Mm to be two to three steps at their At = 96 min magnetogram 
cadence. For similar spatial scale and At (64 and 128 min) we find shorter lifetimes, by 
a factor of two or so, although this is also the spatial scale where our lifetime estimates 
become incoherent. For all spatial scales we investigate, the flow lifetime is shorter than 
that predicted from a linear scaling of lifetime with length scale, sometimes dramatically so. 
We note that our estimates of lifetimes are rough approximations. Longer-lived flows (e.g., 
supergranular flows) would remain undetected in our lifetime estimates if their asymptotic 
autocorrelation lies below our 1/e cutoff. The longer lifetimes of larger-scale flows is another 
possible explanation of the correlation between flow lifet im e and magnetic held strengths: 
stronger magnetic fields might couple nearby regions of plasma more strongly (making the 
plasma more rigid), leading to larger spatial scales for hows in more strongly magnetized 
regions; and such hows should persist for longer than smaller-scale hows. 

In Figure [03 we show htted how lifetimes for x and y components of the how (denoted 
+ and x, respectively) as functions of ( v ), the how speed averaged over space and time 
for each combination of tracking interval At, spatial binning, and apodization parameter a. 
Plot symbols are color-coded by spatial scale, defined as the product of spatial binning and 
apodization parameter. Generally, higher average speeds correspond to shorter lifetimes. 
While a range of lifetimes was found at each average speed, there appears to be a rough 
upper limit at a given speed, with the peak htted lifetime r* scaling roughly as {v)~ 2 . 
Noting that the product of a timescale and the square of speed yields a quantity with the 
dimensions of kinematic viscosity (length 2 /time), the mean and standard deviation of r*(u) 2 
are (2.7 ± 1.9) x 10 11 cm 2 sec^ 1 and (2.9 ±2.3) x 10 11 cm 2 sec -1 for the x and y components 
of the how, respectively. The limit of how lifetimes shown by the dashed line in Figure fl6l 
1.25(u)~ 2 , would correspond to a kinematic viscosity of 7.5 x 10 11 cm 2 sec^ 1 . While the 
physical signihcance of this value is unclear (if it is even physical — it might be an artifact 
of our methods), it roughly agrees with values of viscosity near the photosphere reported by 
Nesis et al. (1990). 


21 


5.3. Lifetimes of Curls &c Divergences 

The lifetimes of normal curls and horizontal divergences of the flow fields, defined as 

( 4 ) 

( 5 ) 


n • (V fc x u) 
(Vk • u) , 


respectivel y, are also of int erest: v ortical flows can efficiently inject magnetic energy and he- 
licity (e.g., iLongcope and Welschl 120001). and di verging flows are as sociated with flux e mer- 
gence in both simulations (lAbbett et al.l l2000l) and observations (jWelsch et all 120091) . In 
Figure [HI we have overplotted mean absolute values of curls (divergences), averaged over 
tracked pixels in x — y — t datacubes of our flow maps, in red (blue), computed using finite 
difference approximations with the smallest binning, twice the native 0.16 "pixel size, in the 
denominator. For consistency with our definition of spatial scale used above, the denomi- 
nator should be the product of macropixel size and apodization parameter a. This would, 
however, introduce an additional inverse scaling with the variable plotted on the horizontal 
axis. Average values of curls and divergences at our smallest length scales for short At are 
on the order of a few times 1CA 4 sec^ 1 , similar to values reported by Verma and Denker 
(2011). For fixed At, curls and divergences fall by roughly two decades over the approxi- 
mately two decades in length scales we study, implying curl and divergence magnitudes scale 
approximately as 1/length. Average magnitudes of curls and divergences decreased by about 
a decade over the range of At’s in our study, a faster decrease with At than for speeds. 


In Figure [T71 we show autocorrelations of normal curls (red) and horizontal divergences 
(blue) computed in subregions of the full FOV, superimposed on 50 G and 200 G contours 
of |i?Los|> f° r fl° ws estimated from 2 2 -binned data, with a = 4 macropixels and At = 64 
min. In Figure [TSJ we show the fitted lifetimes for curls (left) and divergences (right) ver- 
sus subregion- and boxcar-averaged |£>los| from NFI, and subregion-averaged \B Z \ and |B| 
from SP, where, as above, lifetimes were determined by one-parameter fits to subregions’ 
autocorrelation functions assuming exponential decay. As with the lifetimes of x and y com- 
ponents of the flow, curl and divergence lifetimes exhibit a dependence on both field strength 
and the median occupancy of above-threshold pixels in each 32 2 subregion of 2 x 2-binned 
full-resolution pixels. The rank-order correlation of held strength |B| with curl (divergence) 
lifetime, however, is a statistically significant 0.59 (0.58) for the 59 subregions with occupan- 
cies of 1000 or more (of a possible 1024), implying the lifetime versus field-strength trend is 
independent of any occupancy-lifetime correlation. We repeat that our estimates of |5 los| 
suffer from non-monotonicity in the Stokes V/I ratio, while our spatial averages of \B Z \ and 
| B | ignore the likelihood of small filling factors for magnetic held structures. Despite these 
limitations, the existence of a trend seems clear, and we find a stronger dependence of the 
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lifetimes of curls on field strength than for lifetimes of divergences. We note that, as with 
the lifetime of flows themselves, the fitted dependence of lifetime of curls and divergences on 
| B | appears weaker than that on either \B~\ or |-£>los|- 

Finally, in Figure [T9l we show the lifetimes of normal curls (red x’s) and horizontal diver- 
gences (blue +’s) versus spatial scale, defined as above as the product of macropixel bin 
size and windowing parameter cr, for each tracking interval At. Because we estimated flows 
for several macropixel sizes and cr’s, we have multiple samples at some resolutions, most of 
which roughly agree. At most cadences, the lifetimes increase with spatial scale, and then 
saturate. In Figure fl9l we also show linear fits to the logarithms of lifetimes as functions of 
logarithms of spatial scales over the range of spatial scales before saturation occurs. 


6. Summary &; Discussion 

We used FLCT to estimate flows from a ~ 13 hour sequence of high-cadence (~ 2 min), 
high-resolution (0.16 "pixel, 0.3 "diffraction limit) linc-of-sight magnetograms of NOAA AR 
10930. This seeing-free, high-resolution, long-duration dataset enabled us to conduct a 
detailed study of the effects of tracking parameter selection on the estimated flows, and to 
investigate the lifetimes of flows and flow structure (normal curls and horizontal divergences) 
as functions of spatial scale and magnetic held strength. For context, we also characterized 
the lifetimes of magnetic structures. 

Several conclusions can be drawn from our results, which may be divided into practical 
conclusions relevant to the process of how estimation from magnetograms, versus physical 
conclusions about surface magnetic hows and evolution. Our practical conclusions are: 

• Flows estimated by magnetic tracking (and related optical how) methods with a par- 
ticular set of tracking parameters (i.e., tracking interval At between images, pixel size, 
apodization window size cr) are primarily sensitive to evolution on a particular spatial 
scale and time scale. 

• If the tracking interval At between pairs of magnetograms is too short, then changes 
in the held between magnetograms is mostly due to noise instead of actual magnetic 
evolution. This is the noise-dominated regime. Its hallmark is rapid fluctuations in 
estimated velocities, which can be quantified by autocorrelating maps of the velocity 
components. The effect of noise can be reduced by: increasing At, averaging successive 
magnetograms prior to tracking, or both. 

• Because magnetic features persist for much longer than the intensity features to which 
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track ing methods have also been applied (e.g., iBerger et al.l Il998l: IVerma and Denker 
2 0111 ). the upper limit for At is not strongly constrained in magnetic tracking. If, 


however, the tracking interval At between pairs of magnetograms is longer than the 
flow lifetime at that spatial scale, then the estimated flows reflect an average velocity 
over more than one flow lifetime. This is the displacement-dominated regime. Although 
this time-averaged velocity is physically meaningful, it does not represent the actual 
plasma velocity held at a particular time (e.g., the center of the tracking interval). 


• Hence, choice of At can determine whether how estimates are noise-dominated (At too 
short), displacement-dominated (At longer than the how lifetime), or representative 
of actual hows. To estimate actual plasma velocities, appropriate choice of At can be 
established through post-facto analysis of how maps, by autocorrelating hows to find 
a At that both maximizes frame-to- frame autocorrelations (which thus minimizes the 
effects of noise) and is less than the how lifetime r, which we take to be the time at 
which the autocorrelation drops to 1/e. 


• Average estimated speeds decrease with increasing tracking interval At between mag- 
netograms (see Figure fl4l). 


Our physical conclusions are: 


As noted by previous researchers (e.g.. IBerger et al.ll998i ). photospheric hows operate 
over a range of spatial scales, and exist for a range of lifetimes. Hence, referring to 
“the how held,” by itself, is imprecise: explicit reference to spatial and temporal scales 
must also be made. 


• Magnetic structure persists substantially longer in stronger-held regions generally, and 
longer in regions of strong vertical held \B Z \ than in regions of strong magnetic held 
strength |B|. 

• Magnetic structure persists substantially longer when averaged over larger spatial 
scales. Flow components and their derivatives estimated on larger spatial scales also 
persist for significantly longer, though lifetimes scale less than linearly with length 
scale. 


• Flow components, horizontal divergences, and normal curls persist somewhat longer in 
regions with stronger helds: our estimates of their lifetimes are significantly statistically 
correlated with spatially averaged values of each of |£>los|, \B z \, and |B|. We find the 
lifetimes of curls exhibit a slightly stronger dependence on held strength than either 
hows components or horizontal divergences. 
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• Flows with faster average speeds typically exhibit shorter lifetimes. Also, there appears 
to be a rough upper limit on flow lifetime at a given average speed, which scales as the 
inverse of the square of the average speed. 


In this study, we found sets of differing velocity fields for a single magnetogram dataset 
by applying a single correlation tracking method over various spatial and temporal scales. 
This is not surprising: coherent flow structures on varying scales are known to be present 
in the photospheric flow field, e.g., global-scale Sun’s differential rotation (in our synodic 
frame) and meridional flow, high-speed granular flows and lower-speed supergranular flows 
driven by convection, and possibly flows arising from Lorentz forces in magnetized regions. 
In addition, turbulent flows are likely present. This, however, raises the question, “Given 
that a family of velocity fields can be derived on several scales, how can one determine the 
best estimate of the velocity field?” But what determines “best?” It is tempting to try 
to estimate the “true” velocity field in each resolution element (pixel), but what does this 
mean physically? Assuming the photospheric plasma behaves exactly as a fluid down to 
some kinetic scale, then within our smallest resolution element there is a distribution of 
subresolution velocities, and a meaningful definition of the “true” velocity must be defined 
in terms of this distribution. If the distribution of velocities is peaked in speed and direction, 
then some statistical property of each (e.g., mean or median) might meaningfully quantify 
the flow. It is, however, probable that sub-resolution flows are consistent with a power-law 
distribution of flow energy with spatial wavenumber, similar to that proposed by Kolmogorov 
(though perhaps with a different wavenumber scaling), through some inertial range down 
to a dissipation s cale. A reasonable expectation based upon o ur results, as well as those of 
previous studies (IBerger et al.l Il998l: IVerma and Denkerl 1201 11) , is that flows on successively 
smaller scales will have higher speeds. In this case, what we infer from the motion of 
magnetic flux (or of brightness variations, if tracking intensity) is an “effective flow” — the 
combined effect of flows acting on all scales. But it is probable that the effective flow does 
not correspond to “the” physical velocity. 


From the practical standpoint of understanding how flows are related to other phenomena 
of interest (e.g., flares, CMEs, coronal heating, the dispersal of active region magnetic flux), 
we must investigate relationships between these phenomena and flows on varying spatial and 
temporal scales to understand which (if any) scales are most relevant. A key challenge is that 
the largest velocities are present at the smallest scales: while in a sense these strongest flows 
are the easiest to detect, they also have the shortest lifetimes, and are therefore probably the 
least relevant for some processes with longer characteristic time scales, e.g., the emergence 
of large-lcngth-scale magnetic flux systems in active regions. Computing the Poynting flux, 
for instance, based upon an instantaneous estimate of flows on granular scales will likely 
give an incomplete (if not misleading) picture of how much and where magnetic energy is 
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being transported across the photospheric layer. Hence, much work remains to understand 
photospheric flows and their relationships to processes like the transport of magnetic energy 
across the photosphere. 


We remark that the techniques we have applied here can also be applied to both real and 
simulated vector magnetogram sequences, to better understand evolution of both magnetic 
fields and flows. Long sequences of vector magnetograms (from, e.g., HMI) could be tracked 
to analyze flow properties over a range of spatial and temporal scales as functions of the 
vector magnetic field’s properties^ The lifetimes of magnetic and flow fields from simulations 
of photospheric evolution (e.g.. iCheung et ah 2 010 ) can also by analyzed via autocorrelation 
for comparison with lifetimes inferred from observations. 

We e xpect observ a tions using the ne west generation of gr ound-bas ed solar telescop es 


NST Goode et al. 


j201C)ll. GREGOR (Volkmer et al. 2010 1. ATST ( Rimmele et al. 2010 1. 

and EST (Izuccarello and Zuccarello 2011 1 — to play a key role in future studies of flows at 


high resolution. 
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A. Calibration by Co-alignment with SP Vector Magnetograms 


Schriiver et ahl (120081 ) used SP data to produce vector magnetograms of AR 10930 over 12 - 
13 Dec. 2006, which are available online. One of the magnetograms they produced overlaps 
with our tracking run. Here, we describe the procedure we used to interpolate magnetic field 
data from that magnetogram onto a downsampled, cropped grid corresponding to our NFI 
observations. Note that our ultimate goal is to be able to relate the lifetimes of magnetic field 
structures and flows over regions of the NFI FOV to average magnetic properties in those 
regions. Hence, we do not seek pixel-scale agreement between the SP and NFI magnetograms. 
Rather, we seek approximate agreement on the scales of a few arcseconds. 
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Schriiver et al.l (12008 ) report sampling the SP data into 0.63 "pixels. We resampled NFI 


observations of Blos at the start, middle, and end of the 45-minute SP scanning interval from 
the original 0.16 "pixels into (4 x 4) bins to achieve commensurate resolution. Comparison of 
features in the B z map from SP (Figure [201 upper left) with binned Blos maps from NFI at 
the start and end of the scan (Figure [20] middle and bottom of left column) shows that most 
of the discrepancies between the SP and NFI data are in the horizontal direction: features 
in the SP image extend over N Xt sp = 295 pixels in the horizontal direction, while the NFI 
images of the corresponding features are only N Xj nfi = 240 pixels in horizontal extent. 


Given the position of the AR near 20° W from disk center at the time of the scan, and the 
~ 10° extent of the AR, the cumulative effect of horizontal distortion in pixel lengths from 
foreshortening in NFI pixels compared to reprojected SP pixels over the FOV should account 
for a ~ 9% difference, or about 22 pixels. Hence, additional distortion is present. We also 
note that, even if reprojection fully accounted for the SP/NFI discrepancies, un-reprojecting 
the SP data would require a difficult inversion of the convolved effects of solar rotation, 
spacecraft pointing, and spectrograph rastering over the 45 minute scanning period. 


We therefore adopt an ad hoc approach, using repeated Fourier interpolations to spatially 
interpolate the SP data onto each column of the NFI grid. We first interpolated the SP 
data onto N x> nfi evenly-spaced shifts N Xj sp/N Xj nfi apart (one Fourier interpolation per 
shift), but found systematic residual horizontal distortions were present, consistent with a 
higher-order variation in distortion across the image. We therefore included a quadratic 
term in the spacing between interpolations: where Xi is the location of the Ath column of 
NFI data, we interpolated the SP data at x\ = N X) sp/N Xt MFi{ciXi ± hxf). By trial and error, 
we found a = 0.925 and b = 0.06 gave good qualitative agreement between ±100 G and 
±400 G contours of Blos from NFI at the middle of the scan interval with B z from the SP 
magnetogram, as shown in the upper right panel of Figure [201 It may be seen that the NFI’s 
contours of Blos d° not match SP’s B z everywhere, but the remaining offsets do not appear 
to be systematic across the image, or along columns. (We note that the non- monotonicity in 
the NFI signal precludes use of some quantitative measures of NFI-SP agreement [e.g., cross- 
correlation, or minimized square difference] since optimal alignment probably corresponds 
to the weakest umbral held in NFI being co-aligned with the strongest umbral held in SP, 
while features outside of the umbra should generally correspond one-to-one.) 

As discussed above, we are primarily interested in average magnetic properties (vertical 
held strength, total held strength) in larger-scale regions of the NFI FOV, so we only use the 
interpolated SP data binned into (8 x 8) and (16 x 16) macropixels. Comparisons of (32 x 32)- 
binned NFI |Blos| data with (8 x 8)-binned |B.| (and |B|) from SP in the middle-right (and 
bottom-right) panel of Figure [20] shows good spatial agreement at that resolution. 


Note that the SP data do not fully cover the NFI FOV, which is 256 2 pixels when binned 
(4x4). Consequently, we can only use the SP data on a subset of the full NFI FOV. 
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Fig. 4. — The grayscale image is |Blos| in a full-resolution magnetogram, with saturation at ±4 kG. Red 
curves show rank-order autocorrelations of BloSj in (32 x 32)-binned subregions of (2x2)- binned data, 
vs. lag time, up to the duration of the dataset. The vertical range in each cell is [-0.5, 1.0], with a dashed 
black line at zero correlation. Blue curves show one-parameter fits to the decorrelation in each subregion, 
assuming exponential decay, with the decay constant as the only free parameter. Only subregions with 
median occupancy of at least 20% of pixels above our 15 G threshold were fit. As expected, field structures 
persist longer in stronger-held regions. 
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Fig. 5. — Fitted lifetimes of boxcar-averaged |Blos| a nd |B|, averaged over each subregion, for the non- 
bubble region of the FOV. Blue triangles correspond to (32 x 32)-binned subregions (from fits in Figure EP of 
|.Blos|- Bed +’s correspond to (16x 16) -binned subregions |Blos| (autocorrelations and fits at this resolution 
are not plotted in this paper), and green x ’s correspond to |B| averaged over (16 x 16)2 -binned subregions. 
Fitted exponents to log-log values of |Blos| a nd |B| for (16 x 16) bins are shown with corresponding colors. 
Dashed lines correspond to uncertainties in each fitted slope, based upon x 2 : the slope of the blue fit (0.68) 
is similar to the red fit (0.69). 


- 34 - 


8 


CD 

O 

O 


O 

o 


0 

"0 

o 


2*: 

0 

D 

cr 


Frame — to — Frame Corr. w/Varying At (A 2 bin, u = 4 pix.) 

i o F 1 T 1 T ' 1 

Raw = dashed; Avg. = solid 

2 min. = black; 4 min. = aqua; 8 min. = blue; 16 min. = red 

0.8 - 
0.6 \~ 

0.4 

0,2 

0.0 

— 0.2 i i i i i i i i i i i i i i i i i 

0 


Raw = dashed; Avg. = solid 

2 min. = black; 4 min. = aqua; 8 min. = blue; 16 min. = red 



2 4 6 8 1 

Initial Frame Number 


Fig. 6. - Plotted lines show frame-to- frame correlation coefficients for the x- and y- components of estimated 
flows , for varying time intervals At between tracked magnetograms — 2 min in black , 4 min in aqua, 8 min 
in blue, 16 min in red. Dashed curves are from flows derived from un-averaged magnetograms, while solid 
curves are from flows derived from boxcar-averaged magnetograms. Flow estimates are more stable with 
averaged magnetograms and with longer At, until saturation occurs around At = 8 min. 
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Fig. 7. - Successive How maps in a subregion of the full Held of view, derived from unaveragecl magnetograms 
(binned 2x2) separated in time by At = 2 min with o = 4 pixels. Initial Hows are shown by red vectors, 
subsequent Hows are shown by blue vectors. Note signihcant variations between the How Helds in this 
noise-dominated regime. 
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Fig. 8. — Successive flow maps from the same subregion, derived from averaged magnetograms (binned 2x2) 
separated in time by A t = 8 min with a = 4 pixels. Initial flows are shown by red vectors, subsequent flows 
are shown by blue vectors. Note significant agreement between flow fields. 
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Autocorrelations vs. At (2 2 bln) 



Fig. 9. — Plotted lines show correlation coefficients for the x- and y- components of flows estimated at U 
and ti + AT, for several time intervals At between tracked magnetograms — 4 min in black, 8 min in aqua, 
16 min in blue, 64 min in red. Dashed curves are from flows derived with cr = 4 pixels, while solid curves 
are from flows with o = 8 pixels. All flows were derived from boxcar-averaged magnetograms. Horizontal 
dotted line is 1/e. 
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Fig. 10. - Top: Autocorrelations of x- and y-components of flows derived using two values of a (dashed = 
4 pix., solid = 8 pix.) from magnetograms binned into different-sized macropixels (black = 2 2 , blue = 4 2 , 
red = 8 2 ). Evidently, keeping o fixed but doubling bin size approximates keeping bin size fixed but doubling 
a. Additional binning did not improve correlations. Bottom: Autocorrelations of x- and y-components 
of flows estimated with large values of a (16 pix. for all flows here) and from data rebinned into larger 
macropixel sizes (colors as in top panel) for At = 128 min (solid) and At = 256 min (dashed) are long-lived. 
As in top panel, additional binning did not improve correlations. All flows in both panels were derived from 
boxcar-averaged magnetograms. 
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Fig. 11. - The grayscale image is |i?LOs| in a full-resolution magnetogram, with saturation at ±4 kG. 
Autocorrelations of u x (red) and u y (blue) in subregions of the full field of view. The flows were derived 
with At = 64 min and o = 4 pix with (2 x 2) binned data. The vertical range in each cell is [-0.5, 1.0], 
with a dashed black line at zero correlation. Autocorrelations are essentially random in weak-field regions, 
in which few pixels were above the threshold for tracking. 
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Fig. 12. — Fitted How lifetimes for non-bubble subregions of Figure fTTI in which at least 20% of the 32 2 
(2x2) pixels exceeded the tracking threshold, as a function of the spatial average of the unsigned, boxcar- 
averaged LOS magnetic Held in each subregion. Lifetimes for u x in subregions are shown with red +’s, u y 
are blue x’s, and Hts to these are red and blue solid lines. The red A’s (green V's) are lifetimes of u x 
versus subregion- aver aged \B Z \ (|B|) from the SP data, and the black (green) line is a Ht, with slope 0.17 
(0.13). The correlation between How lifetime and Held strength is statistically signiHcant, though the Htted 
power-law dependence is weak. 
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Fig. 13. — Black numerals show average flow speed (left axis) vs. spatial scale , parametrized by the product of 
the macropixel bin size and the windowing parameter cr, for several flow cadences; the numbers correspond 
to log 2 (bin size), in units of 0.16 "pixels. Red (blue) numerals show the log of the average normal curl 
(horizontal divergence) of the flow, on the right scale, and are slightly offset right (left) for clarity. Because 
we estimated flows for several macropixel sizes and cr’s, we have multiple samples at most resolutions. For 
shorter At, bin size affects the average flow speed at a fixed spatial scale; for longer At, bin size does not 
strongly affect average flow speed at a fixed spatial scale. In contrast, binning seems to affect average flow 
curls and divergences more at larger At. 
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Avg. Speed at each At 



Fig. 14. — Flow speed as a function of cadence, averaged over both space and time for each spatial resolution 
(binning xcr) and over spatial resolutions at each At. Error bars show standard deviation over resolutions. 
Least squares linear fit to the logarithms of cadence and speeds returns a power-law exponent of -0.34. 
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Fig. 15. — Flow lifetimes (+’s for u x , x’s for u y ) as a function of flow spatial scale, as parametrized by 
the product of the macropixel bin size and the windowing parameter o, for several flow cadences. Because 
we estimated flows for several macropixel sizes and cr’s, we have multiple samples at some resolutions. We 
computed linear fits (solid lines) to the logarithms of spatial scale and lifetime over a limited range of scales. 
The fitted power-law exponent is printed in the title of each cadence’s plot. 
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Fig. 16. — Fitted lifetime vs. flow speed averaged over space and time for each choice of At, spatial binning, 
and o. Lifetimes for u x and u y are plotted with + ’s and x ’s, respectively. Generally, higher average speeds 
correspond to shorter lifetimes. While a range of lifetimes exists at each average speed, there appears to be 
an upper limit at a given average speed, with the peak fitted lifetime scaling roughly as ( average speed)~ 2 
(dashed line; note that this line is not a fit). Points are color-coded by spatial scale of the flow (binning xo). 
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Fig. 17. - The grayscale image is |7?los| in a full-resolution magnetogram, with saturation at ±4 kG. 
Autocorrelations of normal curls (red) and horizontal divergences (blue) in each subregion of the full field 
of view are overplotted. The flows were derived with At = 64 min and a = 4 pix. The vertical range in 
each cell is [-0.5, 1.0], with a dashed black line at zero correlation. Autocorrelations are essentially random 
in weak-field regions, in which few pixels were above the threshold for tracking. 
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Fig. 18. — Fitted lifetimes of normal curls (left) and horizontal divergences (right) for subregions of Figure 
E2 in which at least 20% of the 32 2 (2 x 2) pixels exceeded the tracking threshold, as functions of the 
spatial average |!?los| (from NFI; blue +’s), \B Z \ (from SP; red A’s), or |B| (from SP; green V’sj in each 
subregion. The correlation between curl/divergence lifetime and field strength is statistically significant, 
though the fitted power-law dependence is weak. The magnetic dependence of the lifetimes of curls is 
somewhat stronger than that of lifetimes of divergences. 
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Fig. 19.— Lifetimes of normal curls (red x’s) and horizontal divergences (blue +’s) as a function of flow 
spatial scale, as parametrized by the product of the macropixel bin size and the windowing parameter o, 
for several flow cadences. Because we estimated Hows for several macropixel sizes and o ’s, we have multiple 
samples at some resolutions. We computed linear fits (solid lines) to the logarithms of spatial scale and 
lifetime over a limited range of scales, and list the Htted power-law exponents on each cadence’s plot. 
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Cropped SP, : 2006 Dec. 12, 20:30 - 21:15 Grayscale: Interpolated SP; Countours: NFI 
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Fig. 20. - Left column: a) SP map of B z , scanned over 45 minutes, b) & c) NFI images binned to (4 x 4) 
from 0.16 "pixels, from the beginning and end of the SP scan. Note lateral extent of SP data is 295 pixels, 
while NFI data are 240 pixels. Hence, the primary difference between the datasets is lateral stretching. Right 
column: d) Contours of Blos from scan-centered NFI magnetogram (—/+ 100G and 400 G in red/blue, 
resp.) overlain on the interpolated SP data. Slight inaccuracies in co-alignment vary across the FOV. e) & 
f) |i?LOS I from NFI contoured over grayscale of\B z \ and |B| from SP, resp., both binned (8 x 8). Contours 
of | B LOS I from NFI correspond to 100, 200, 400, 800, and 1600 G. 








